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Abstract. 

The extended Bose-Hubbard model subjected to a disordered potential is predicted 
to display a rich phase diagram. In the case of uniform random disorder one finds 
two insulating quantum phases - the Mott-insulator and the Haldane insulator - in 
addition to a superfiuid and a Bose glass phase. In the case of a quasiperiodic potential 
further phases are found, eg the incommensurate density wave, adiabatically connected 
to the Haldane insulator. For the case of weak random disorder we determine the 
phase boundaries using a perturbative bosonization approach. We then calculate the 
entanglement spectrum for both types of disorder, showing that it provides a good 
indication of the various phases. 
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1. Introduction 

The extraordinary experimental advances on the realization and control of ultracold 
quantum gases subjected to optical lattice potentials |1] pave the way to the application 
of these systems as 'quantum simulators' |2], capable of exploring with unprecedented 
accuracy complex models from condensed matter physics [Sj H] to high-energy physics 

The investigation of the interplay of disorder and interaction effects remains one 
open question in condensed matter physics, linked to the study of the metal-insulator 
transition. For bosonic systems, the Bose glass phase [HI [7] is an example of a novel 
strongly correlated phase arising from the simultanous effect of disorder and interactions. 
Most of solid-state based physical systems have to deal with some amount of disorder, 
which originates e.g. from defects in the material or impurities atoms. A very peculiar 
feature of quantum gases is that one can add a tunable and controllable amount 
of disorder in the pure system. In the regime of vanishing interactions, Anderson 
localization has been observed, first in one spatial dimension [HI |9] and then also in 
three dimensions [TOl [TT] . The Bose glass phase has also been explored with bosons 
with short-range interactions on a lattice [12] . 

The very recent advances on trapping and cooling ultracold molecules [13] and 
atoms with a large dipole moment [TH [15] prelude to the exploration with atomic 
quantum simulators of yet another set of systems, those with long range interactions. 
For the case of one-dimensional bosons on a lattice, the minimal model accounting for 
longer-range interactions is the extended Bose-Hubbard model (EBHM), which includes 
next-neighbour interactions among the bosons. In the case of a clean system, the EBHM 
displays already a rich phase diagram, which features novel several insulating phases, ie 
the density wave and the Haldane insulator, in addition to the Mott insulator already 
found in the Bose-Hubbard model [7]. 

In this context, a relevant question arises on the fate of such insulating phases under 
the effect of disorder, as well on which novel phases arise in the presence of disorder. 
In the strong-coupling regime and for unitary lattice filling, a good starting point to 
understand this behaviour is to map the problem onto spin chains. For the latter system, 
the stability of the various insulating phases at weak disorder has been studied using 
renormalization group arguments [16]. For example, the density wave phase is found to 
disappear at infinitesimally small disorder, according to the Imry-Ma argument [TTlfTS]. 
In analogy to the results known for the Bose-Hubbard model, other insulating phases are 
expected to shrink in the phase diagram, in favour of disordered correlated phases of the 
Bose-glass type. The mapping onto spin systems fails when the tunnel energy becomes 
sufficiently important with respect to repulsive interactions to allow large occupancy of 
given lattice sites, and a superfiuid, gapless phase builds up. In the general case then 
the problem has to be addressed numerically, and we have recently explored it using 
Density-Matrix Renormalization Group [19], establishing a phase diagram by following 
the behaviour of several correlation functions. In this work, we complement and support 
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the phase diagram using two complementary tools: a bosonization and renormalization 
group description at weak disorder, and the study of the entanglement spectrum for the 
system. 

2. The model, numerical methods and phase diagrams 

We consider a system of N dipolar bosons confined onto a one-dimensional deep optical 
lattice and in the presence of a very shallow trapping potential. We assume the 
dipole orientation to be perpendicular to the lattice direction and truncate the dipole- 
dipole interaction potential to nearest-neighbour interactions. Although longer range 
interactions do play a role, for very weak interactions and sufficiently large dipoles, the 
most relevant properties of polar physics can be already understood from a model of 
nearest neighbour interaction. This leads to the Hamiltonian of the extended Bose- 
Hubbard model (EBHM) 

U ^ 

H = -t^{blbi+i + h.c.) + — ^rii^Ui - 1) + V^nini+i + ^e^ni, (1) 

i i=l i i 

where b\, bi are creation and annihilation operators for bosons at site i, ui = b\bi are the 
number operators, t is the tunnel energy, and U {V) are the on-site (nearest-neighbour, 
NN) interaction energies. Disorder is included in the model either through the random 
on-site energies e^, chosen to be uniformly distributed in the interval [—A, A], or by a 
quasiperiodic potential = A cos(27rai + (fy and a is a rational number. 

2.1. Numerical method 

We have studied the EBHM using the DMRG method with open boundary conditions. 
The considered system sizes range up to 233 sites and we have taken up to 60 disorder 
realizations per point. In order to avoid the presence of metastable states we allow 
the number of optimal states to shrink or expand at every DMRG step according to 
a two-step algorithm which keeps at least one of the eigenvectors in the blocks of the 
reduced density matrix if they have only zero eigenvalues, and then keeps an additional 
eigenvector with zero eigenvalue in the block with non-zero eigenvalues if the eigenvalues 
decay very sharply to zero. Furthermore, we have eliminated the edge states in the HI 
phase by adding one more particle or by coupling two extra hard-core bosons at the 
edges of the chain in order to form a singlet state [2Dj . 

2.2. Phase diagram of the EBHM without disorder 

The phase diagram of the EBHM in absence of disorder is known [2T1 and is 
illustrated in FigJU If the tunnel energy is dominating on the on-site interactions 
the bosons are delocalized throughout the lattice and the system is superfluid. The 
fluid displays quasi-long range order, with an algebraic decay of the first-order 
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Figure 1. Phase diagram for polar bosons witli nearest-neiglibor interactions. Left 
panel: clean case in the plane (V/t, U/t), with inset showing the value of the Luttinger 
parameter Ka along the MI-HI critical line, which is described by a Luttinger model. 
The HI-DW critical line belongs to the Ising universality class [26]. Central and right 
panel: case of uniform disorder, in the plane {V/t, A/t) for U/t = 5 and U/t = 3. 

correlation function G{r) = (blbi^r) / a/ (n-i) (rij+r) oc r'^^"^^ . At increasing on- 
site interactions and weak NN interactions a Kosterlitz-Thouless transition occurs 
towards the incompressible Mott-insulator (MI) phase, characterized by a hidden 
parity order Op = lim|j_j|^oo((-l)^'<'<^ '^'^0 i22j, , where Sni = 1 - m. At 
sufficiently large values of U and intermediate values of a second insulating phase 
is found. This Haldane insulator (HI) is characterized by a hidden string order 
Os = \im\i_j\_^ao{Sni{—l)^i<'<i^"''^6nj). At sufficiently large values of V and U a 
third insulating phase occurs, a density wave with spatial modulation (DW) as can 
be identified by the finite correlator Ocw = \im\i_j\^oo{{—^y~'^ ^niduj) . At sufficiently 
large U a direct ffist-order transition is expected between Mott-insulator and density- 
wave. 

2.3. Phase diagram of ttie disordered EBHM 

The phase diagram of the EBHM in presence of uniform disorder has been ffist studied 
in [I9] and is illustrated in FigJH Among the phases of the pure system, the density wave 
phase is unstable under infinitesimal disorder according to the Imry-Ma argument [T7] 
thus disappears. Mott insulator and Haldane insulator disappear at sufficiently strong 
disorder to leave place to a compressible, non-superfluid Bose-glass phase [231 El Ej- 
An additional superfluid lobe is found at finite disorder as in found in Ref. [Ml [25] for 
= 0, and corresponds to the regime where repulsions overcome localization effects [6]. 
It is not easy to resolve numerically the behavior of the critical MI-HI point for 
weak disorder, and in particular whether the critical line remains stable or whether a 
Bose glass intermediate region opens at A = 0. This issue is discussed in Sec. [3] using 
bosonization and renormalization (RG) arguments. 

2.4. Phase diagram of the EBHM with quasiperiodic potential 

For the case of a quasiperiodic potential the phase diagram, shown in Figj2l considerably 
differs from the one with uniform disorder [19]. The main features are the presence of 
an incompressible density wave phase, typical of the quasiperiodic potentials [271 EH], 
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adiabatically connected to the Haldane insulator phase, and the persistence of a density 
wave phase, due to the intrinsecally different nature of the quasiperiodic potential with 
respect to a truly random potential [29] . 



3. Bosonization approach at weak disorder 

For sufficiently strong interactions, where number fluctuations on each site are relatively 
small, we truncate the occupancy of each site to the values {0, 1,2}. We then employ 
the Holstein-Primakoff transformation = dui = 1 — nj, Sf = y/2 — Uihi. to map the 
EBHM onto a spin-1 Hamiltonian with single- ion anisotropy 

H=-2tY, S:S!^, + SfSl, + VJ2 S^SUi + U/2 J^iSn' (2) 

i i i 

Following the early works of Timonen and Luther and Schulz [30l l3T] we represent the 
spin-1 operators as the sum of two spin 1/2 operators, S"" = af ^ + af2- This brings the 
Hamiltonian ([2]) into the one of two coupled spin chains. Furthermore, a Jordan- Wigner 
transformation is employed to map the spin-1/2 operators onto fermions according to 
'^ix — x^i>' ~ 1/2) '^tx ~ x^^'^'^^^^^li x^ri,x, with A = 1,2. A continuum limit 
0'n,x = V(^J2p=±'^pi^^) where a is the lattice spacing and p stands for ±, left or right 
mover is taken. Finally, one employs a low-energy description of each fermionic field, 
ipp{x) ~ ^i-e^J'^-p'^e"**^^'^^*^^''"^^*^^)'', where the fields 6\{x) and (px^x) satisfy canonical 
conjugation relations [(I)\{x),dx6\i{x')\ = m5\\'5{x — x'). We have kp = n/{2a) when 
(S^) = i. e. the filling is one boson per site. This leads to the Hamiltonionan of two 
coupled Tomonaga-Luttinger fluids [311 [Ml [IE] , which takes a simple form H = Ha + Ho 
once the 'acoustical' and 'optical' combinations are introduced (pa = {<Pi + 02)/\/2, 
and similarly for the 9x fieds. 



Ha = ^ [ dx[Ka{dJaf + -^{dM'] + 7^ / dxcos{V8<Pa) (3) 
2tt J Ka [Tray J 

Ho = ^ [ dx[Ko{d^eof + ^{dAof] + f dxcosiVs^Po) 

2tt J Ko (vra)^ J 



+ 7^ / dxcos{V2eo) (4) 



where gi = g2 = {U — V)a, gs = —2ta. Coupling between acoustical and optical sectors 
is found at higher order [2T| [22] and is therefore less relevant than the terms listed here. 
The weak-coupling expressions for the Luttinger parameters entering Eqs.(|3]) and 
read 



Ua = 2ta\l 1 + Ka = . (5) 

2nta 1^ U+&V 



2iTta 



Uo = 2ta\l 1 - ^ ^ , Ka = (6) 
27Tta /i _ u-v 

2nta 
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Introducing the rescaled fields 0+/- = 9a/o/V^ and = \/20a/o, the Hamiltonians 
in Eqs. ([3]) and (jl]) can be brought to the form used in Refs. [211 122] • 

The phase diagram [31] can be deduced from ([3]) and (jlj). For Ka > 1 in (|3]), the 
cosine term is irrelevant, and the 'acoustic' modes are gapless. For Ka < 1 the cosine 
term is relevant, and the field (j)a is pinned either to for < or to tt/a/S for gi > 0, 
the spectrum of Ha being always gapped. Meanwhile, in @, at least one of the two 
cosines is relevant so the spectrum is always gapful. Depending on the parameters, 
either 6^ or (po is pinned. When 0o (resp. 6o) is pinned, the correlation functions of its 
dual fields e*^^° (resp. e*^"^") decay exponentially. Combining the different possibilities, 
we obtain the two gapless superfluid phases (when Ka > I and either 0o or 9o is pinned), 
the gapful density wave phase (when both and (po are pinned), and two phases with 
(pa pinned and 9o pinned. The phase with (pa pinned to zero (i. e. gi < 0) is the Mott 
insulator, while the phase with (pa pinned to 7r/\/8 (i. e. gi > 0) is the Haldane insulator 
[2I1I22]. 

We focus in the following on the parameter regime corresponding to the Haldane- 
insulator to Mott-insulator transition point at weak disorder. This regime is difficult to 
access numerically, but is amenable to a perturbative renormalization group calculation. 
At the critical point gi = which separates the Mott-insulating from the Haldane 
insulating phase, the system is in a Luttinger liquid state [321 [331 EH [35] and the boson 
Green's function decays as 

1 

(6t(x)6(0)) = \C\' (^^j . (7) 

One kind of randomness that we are able to treat is the on-site disorder of the form 
Yln^nbnbn, which for the spin-1 chain corresponds to the effect of a random field along 
the z-axis [361 [16]. When dealing with the coupling to disorder, one has to bear in mind 
that the bosonized expression of the boson field is a series which contains higher order 
harmonics |3 7] in the field (pi 2 from which the boson number operator can be written as 



b\firi = Am cos{mV2(pa — 2mkFx) cos(m\/20o), (8) 

m=l 

where x = na. The term of order one has been treated in [16] and it is relevant when 
Ka + Ko < 3, while higher order terms have been neglected. When the second order term 
is taken into account, A2 cos{2 \/2(p a — 4:kFx) cos{2y/2(po) , and performing a perturbation 
theory in cos{2y/2(po) at first order, an effective coupling to disorder is generated of the 
form: 

H:ff = (W(a:)A2e*^'^'' + h.c), (9) 



where ^4^^,(2;) is a gaussian random variable with ^4:kp{x)^li^p{x') = Deff6{x — x'). A 
renormalization group treatment of such a term gives 



dl 



i3-AKa)D,ff, (10) 
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Figure 2. Phase diagram for polar bosons witli nearest- neiglrbor interactions in a 
quasiperiodic potential in the plane {V/t, A/t) for U/t = 5. Left panel: phase diagram 
from the study of correlations functions as in |19j . Right panel: results from the 
entanglement spectrum. 

and Deff is relevant when Ka < 3/4. The terms with m > 2 become relevant for lower 
values of Ka < Thus, along the critical line = one expects a stable Luttinger 

liquid line till = 3/4, and below it a Bose glass phase between the Mott-insulating 
and Haldane-insulating phase takes place. 

Let us note that even in a model where only the m = 1 term is kept in the expansion 
dH]), the term ([9]) will be generated [36] by integrating over the fluctuations of the field 
9o- So the condition of stability K* > 3/4 is a generic one. As a consequence of this 
analysis, we expect that the HI-MI transition point in the presence of disorder has a 
'Y' shape for Ka > 3/4 and a 'V shape for Ka < 3/4. This is in agreement with the 
numerical results shown in FiglH where for U/t = 5 with Ka = 0.6 we have a 'V case 
and for U/t = 3 with Ka = 0.8 we have a 'Y' case. 



4. Entanglement spectra 

The calculation of the entanglement spectrum is a novel approach to identify the 
quantum phases in several models. It is particularly useful to bring insight into 
topological phases and phases with nonlocal order parameter (see eg [3H1 EHl SOl 1^ ) 
and has been specifically analyzed for the case of the Bose- Hubbard model [201 1121 112] . 
The entanglement spectrum is defined as the spectrum — InAj(L^) of the effective 
Hamiltonian — Inp^, obtained by partitioning the system density matrix into two parts 
A and B (of length L^+Lb = L) and tracing over the B part. It has been shown that the 
behaviour of the eigenvalues Aj(L^) and their degeneracy differs in the various phases 
of the clean EBHM [20], thus allowing to infer the structure of the phase diagram. We 
show here how the study of the entanglement spectrum can also give useful information 
in the disordered and aperiodic case. 

In order to obtain the phase diagram, we take the combination of the first four 
largest eigenvalues C = Af — + — Aj , where Af = {1/L) Xi{LA)- The result 
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Figure 3. Entanglement spectrum of the EBHM under quasiperiodic potential: largest 
eigenvalues as a function of the partition along the chain corresponding to a) MI at 
V/t = 0.5 A/t = 0.2; b) HI, at V/t = 3.3, A/t = 0.2; c) DW, at V/t = 3.6, A/t = 0.2; 
d) ICDW, at V/t = 1 A/t = 6; e) SF, at V/t = 0.5, A/t = 3; f) BG, at V/t = 2, 
A/t = 7.8. 

is illustrated in Figj2j The alternating regions of low and large values of ( have a very 
good correspondence to the phases predicted from the study of correlation functions. 
This is readily explained from the study of the behaviour of the largest eigenvalues 
Aj(Lyi) in the various phases, of which some examples are given in FigJSl In the Mott 
insulator phase in the clean case a degeneracy is found between and Xj, as well 
as between and A^ while Af is non degenerate. This feature is found to persist 
also in the quasiperiodic case, and yields a large value for (. In the superfluid phase 
instead, the largest eigenvalues are almost equidistant; hence yielding a small (. In the 
clean case, the Haldane insulator phase is characterized by a double degeneracy of the 
largest eigenvalues, hence implying another region of vanishing ( 0. Such degeneracy is 
gradually broken at increasing the strength of the quasiperiodic potential. In the density 
wave phase no degeneracy is found in the clean case, thus allowing to clearly distinguish 
this phase from the neighbouring Haldane insulator; in the presence of a quasiperiodic 
potential, also this phase gradually disappears. A similar level structure is found for 
the incommensurate density wave. Finally, the BG phase displays a structure close to 
the ICDW phase, with a few large non-degenerate eigenvalues. 

We consider next the case of a uniform disorder. In Fig. H] we show the phase 
diagram obtained from the entanglement spectrum. Some examples of the largest 
eigenvalues in the various phases are given in Fig. |5] for a given disorder realization. 
The main features are the same as in the quasiperiodic case: a unique large eigenvalue 
for the MI case and a few equally spaced eigenvalues for the SF case, the BG being 

I Since the simulation in Figl2][b) is done at integer filling the edge states cannot be eliminated, hence 
a small breaking of the double degeneracy is found at vanishing disorder, as in [19] . 
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V/t V/t 

Figure 4. Phase diagram of the uniform disorder case from the entanglement 

spectrum, (a) at U/t = 5, (b) at U/t = 3. Calculations performed with L — 55 
and averaging over 20 disorder realizations. 




Figure 5. Entanglement spectrum of the EBHM in a single realization of the uniform 
random disorder: largest eigenvalues as a function of the partition along the chain 
corresponding to a) MI at V/t = 0.5 A/t = 0.2; b) HI at V/t = 3.3, A/t = 0.2; c) 
BG (remnant of DW) at V/t = 3.6, A/t = 0.2; d) SF V/t = 0, A/t = 3.2; e) EG at 
y/i = 3.1, A/i = 2. The onsite interaction is U/t = 5 and L = 89. 

intermediate between the two above configurations. We notice that although the DW 
disappears under finite-size scahng, some remnants of this phase are still visible in the 
eigenvalues at very low disorder strength. 

5. Conclusions and perspectives 

In conclusion, we have studied the phase diagram of the extended Bose-Hubbard model 
in the presence of a quasiperiodic potential or a uniform random disorder. The analysis 
of the entanglement spectrum is found to be very helpful in identifying the various 
phases obtained from the analysis of correlation functions, and well matches with the 
predictions of a renormalization group approach, when applicable. In perspective, it 
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would be interesting to further explore the nature of the various phase boundaries, 
extending the work done in the case of the Bose-Hubbard model H5] . 

6. Acknowledgements 

We thank T. Vekua for useful discussions. X. D. and L. S. are supported by the German 
Research Foundation (SA1031/6), the German-Israeli Foundation, and the Cluster of 
Excellence QUEST. E. O. and A. M. acknowledge support from the CNRS PEPS-PTI 
project "Strong correlations and disorder in ultracold quantum gases", and A. M. from 
the Handy- Q ERG project. 

References 

[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008). 

[2] Imnianuel Bloch, Jean Dalibard, Sylvain Nascimbene, Nat. Phys. 8, 267 (2012). 

[3] U. Schneider et al., Science 322,1520 (2008). 

[4] R. Jordens et al.. Nature 455, 204 (2008). 

[5] M. Endres et al. Nature 487 454, (2012). 

[6] T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 325 (1988). 

[7] M. P. A. Fisher et al, Phys. Rev. B 40, 546 (1989). 

[8] G. Roati et al, Nature (London) 453, 895 (2008). 

[9] J. Billy et al. Nature (London) 453, 891 (2008). 
[10] S. S. Kondov, W. R. McGehee, J. J. Zirbel, B. DeMarco, Science 334, 66 (2011). 
[11] F. Jendrzejewski et al, Nat. Phys. 8 398 (2012). 
[12] L. Fallani et al, Phys. Rev. Lett. 98, 130404 (2007). 
[13] S. Ospelkaus et al., Science 327, 853 (2010). 
[14] A. Griesmaier et al., Phys. Rev. Lett. 94, 160401 (2005). 
[15] M. Lu et al, Phys. Rev. Lett. 107, 190401 (2011). 
[16] V. Brunei and Th. Jolicoeur, Phys. Rev. B 58, 8481 (1998). 
[17] Y. Imry and S. Ma, Phys. Rev. Lett. 35, 1399 (1975). 
[18] R. Shankar Int. J. Mod. Phys. B 4, 2371 (1990). 

[19] X. Deng, R. Citro, E. Orignac, A. Minguzzi and L. Santos, larXiv:1203.0505l (2012) 
[20] X. Deng and L. Santos, Phys. Rev. B 84, 085138 (2011). 

[21] E. G. DaUa Torre, E. Berg, and E. Ahman, Phys. Rev. Lett. 97, 260401 (2006). 

[22] E. Berg et al, Phys. Rev. B 77, 245119 (2008). 

[23] T. Giamarchi and H. J. Schulz, Europys. Lett. 3, 1287 (1987). 

[24] S. Rapsch, U. Schollwock, and W. Zwerger, Europhys. Lett. 46, 559 (1999). 

[25] N. V. Prokof'ev and B. V. Svistunov, Phys. Rev. Lett. 80, 4355 (1998) 

[26] C. Degh Esposti Boschi, E. Ercolessi, F. Ortolani, M. Roncagha, Eur. Phys. J. B 35, 465 (2003). 
[27] T. Roscilde, Phys. Rev. A 76 063605(2008). 

[28] G. Roux, T. Barthel, I. P. McCuUoch, C. KoUath, U. SchoUwck, and T. Giamarchi, Phys. Rev. A 

78, 023628(2008). 
[29] M. Albert and P. Leboeuf, Phys. Rev. A 81, 013614 (2010). 
[30] J. Timonen and A. Luther, J. Phys. C: Sohd State Phys. 18, 1439 (1985) 
[31] H. J. Schulz, Phys. Rev. B 34, 6372 (1986). 

[32] W. Chen, K. Hida and B. C. Sanctuary, Phys. Rev. B 67, 104401 (2003). 

[33] C. Degh Esposti Boschi, E. Ercolessi, F. Ortolani and M. Roncagha Eur. Phys. J. B 35, 465 (2003). 
[34] A. F. Albuquerque, C. J. Hamer and J. Oitmaa Phys. Rev. B 79, 054412 (2009). 
[35] S. Hu, B. Normand, X. Wang, and L. Yu Phys. Rev. B 84, 220402(R) (2011). 



Bosonization and entanglement spectrum for ID disordered polar bosons 



11 



[36] E. Orignac and T. Giamarchi, Phys. Rev. B 57, 5812 (1998). 

[37] F.D.M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

[38] H. Li and F.D.M. Haldane, Phys. Rev. Lett. 101, 010504 (2008). 

[39] P. Calabrese and A. Lefevre, Phys. Rev. A 78 032329 (2008). 

[40] N. Regnault, B. A. Bernevig, F. D. M. Haldane, Phys. Rev. Lett. 103, 016801 (2009) 

[41] A.M. Lauchli, E.J. Berghohz, J. Suorsa and M. Hacque, Phys. Rev. Lett. 104 156404 (2010). 

[42] V. Alba, M. Haque, A.M. Laeuchh, Phys. Rev. Lett. 108, 227201 (2012). 

[43] V. Alba, M. Haque, A.M. Laeuchh, larXiv:1212.563l 

[44] E. Ahman, Y. Kafri, A. Polkovnikov, and G. Refael, Phys. Rev. B 81, 174528 (2010). 

[45] Zoran Ristivojevic, Aleksandra Petkovic, Pierre Le Doussal, and Thierry Giamarchi, Phys. Rev. 
Lett. 109, 026402 (2012) 



